Electrophoretic separation of proteins via complexation with a polyelectrolyte 
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We suggest to augment standard isoelectronic focusing for separation of proteins in a gradient 
of pH by a similar focusing in the presence of a strongly charged polyelectrolyte (PE). Proteins 
which have the same isoelectric point but different "hidden" charge of both signs in pi point make 
complexes with PE, which focus in different pH. This is a result of charge inversion of such proteins 
by adsorbed PE molecules, which is sensitive to the hidden charge. Hence the hidden charge is a 
new separation parameter. 



Separation of proteins in a gradient of pH, or isoelec- 
tric focusing is the basis of proteomics (|,0]. It uses the 
fact that proteins have both basic and acidic groups. Let 
us assume that on the left side of electrophoretic cell pH 
is so small that all basic groups are positively charged 
and all acidic ones are neutral, so that the protein is net 
positive and moves in electric field to the right where pH 
increases. While protein moves to the right more and 
more basic groups are neutralized and some acidic one 
dissociate and become negative. As a result for each pro- 
tein there is the so called isoelectric point pH = pi, where 
its charge changes sign. Beyond this point a protein be- 
comes negative and returns back. Thus, each protein 
accumulates near its pi. It is not, however, totally sat- 
isfactory because there are many different proteins with 
the same pi. One still needs another parameter to sepa- 
rate them. 

A standard method to do this consists in addition of a 
detergent SDS, which has a long hydrophobic tail and a 
negatively charged head |!|. Molecules of SDS denature 
a protein globule by attaching their tails to hydrophobic 
parts of protein and cover it by charged heads. Result- 
ing rod-like negatively charged complexes are separated 
by their mobility, which in turn depends on the length 
of the protein and the number of its hydrophobic groups. 
This is the second coordinate of two-dimensional elec- 
trophoresis of proteins (pH is the first). 

In many cases even this two-dimensional analysis does 
not provide necessary resolution. Furthermore, native 
structure of proteins is irreversibly destroyed by SDS and 
the protein can not be used for a farther analysis. There- 
fore, if possible, another method should be used together 
with isoelectric focusing. 

In this paper, we suggest a different idea for separa- 
tion of proteins with the same pi. In addition to the total 
charge of protein this method is sensitive mostly to the 
absolute value, q, of the positive and negative "hidden" 
charges at its pi. Our goal is to separate proteins with the 
same pi, but different q. For this purpose, we suggest to 
complex protein globules with a strongly charged, short 
polyelectrolyte (PE). For different proteins with the same 
pi protein-PE complexes have a different number of PE 
molecules which depends on q. Correspondingly, isoelec- 



tric points of these complexes (values of pH where they 
are neutral) differ from a protein to protein. Therefore, 
complexes of different proteins can be separated by iso- 
electric focusing in a gradient of pH. Using short PE is 
important because a long one could bind several different 
proteins together. 

Let us for a given protein plot on Fig. 1 the isoelectric 
points of the protein (thin line) and protein-PE complex 
(thick line) as a function of the concentration of PE, N. 
Adding to this plot the signs of the bare protein and the 
protein-PE complex we obtain the " sign" phase diagram 
in the plane (N, pH) 
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1. The sign phase diagram in the plane (N, pH) for 
the scenarios A (a) and B (b). The horizontal thin line corre- 
sponds to pH = pi, the curve pH = pl(q, N) is shown by the 
thick line. 
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In the case of small enough N, when there are no com- 
plexes at all, plus and minus are separated by the hori- 
zontal thin line pH(iV) = pi. We show below that there 
are two different scenarios for deviation of the thick line 
from the horizontal line, when complexes appear. They 
lead to two different "sign" phase diagrams. According 
to the scenario A (Fig. la), deviations from horizon- 
tal straight line happens when concentration of PE, N, 
reaches the critical concentration N% at which the neu- 
tral protein adsorbs the first molecule of PE. At N > Ni 
this complex is neutral at pH = pl(q, N) > pi. Thus, at 
N = Ni the horizontal line pH = pi, crosses over to pH = 
pl(g, N) curve, which deviates in the direction of larger 
pH. At some larger N — N2 the point ~pl(q, N) reaches 
pH =14. The complex is positive at pH < pl(q,N) and 
negative at pH > pl(<7, N) (Fig. la). 

The scenario B is different from A, because in this case 
Ni > N2 (see Fig. lb). In this scenario again the neu- 
tral protein at pH = pi does not adsorb a PE molecule 
till N = N\. However, at pH = 14 where protein is 
strongly negatively charged the charge of the complex 
can become positive already at smaller concentrations 
N2 < N < N\. As a result at a fixed N in the range of 
concentrations N2 < N < N\ , protein changes sign twice 
with growing pH (Fig. lb). First time this happens in 
the standard pi point of the bare protein. At pH > pi 
protein becomes negative, so that at some point it starts 
adsorbing PE molecules. If pH increases further at some 
point pl(g, N) protein charge becomes neutralized by ad- 
sorbed PE molecules. At even larger pH the bare charge 
of protein is so large that it becomes overcharged by PE. 

In both scenarios the new isoelectric point pl(q, N) de- 
pends both on q and N. In other words, it is different 
for protein-PE complexes of different proteins with the 
same pi. This is the basis for proposed method of separa- 
tion of proteins with the same pi. Although we literally 
defined q as a hidden charge one can also think that this 
notation in formula pl(q,N) includes other hidden pa- 
rameters which discriminate between different proteins 
with the same pi. 

In a simple-minded experiment one can cut out a stripe 
of the gel, where proteins with a given pi and different q 
are focused and put into another device with orthogonal 
to the stripe pH gradient and with the concentration N 
of a PE. In this paper, we have in mind this simple set 
up although it is possible that in future one can create 
a gradient of N orthogonal to the gradient of pH and 
visualize Fig. 1. 

Before explanation of the origin of the scenarios A and 
B let us emphasize the common nontrivial feature of these 
scenarios. Both of them result from the phenomenon of 
charge inversion of a protein by PE. Let us consider it on 
the example of large pH close to 14, when all basic groups 
of the bare protein are neutralized and all acidic ones 
are ionized so that the bare protein has a large negative 
charge. Nevertheless, if the concentration of added PE, 



N, is large enough, the total charge of PE molecules stick- 
ing to the protein can be even larger than bare charge so 
that the net charge of the protein-PE complex is positive. 
In other words, a strongly negatively charged bare pro- 
tein already neutralized by adsorbed PE molecules con- 
tinues to attract new PE molecules. This counterintuitive 
phenomenon is called charge inversion. It becomes pos- 
sible because PE molecules repel each other and form a 
correlated liquid at the surface of the protein. A new PE 
molecule approaching the neutral protein-PE complex re- 
pels already adsorbed molecules on the protein surface 
and creates a correlation hole in this liquid, which plays 
the role of the electrostatic image of PE. As a result the 
new PE molecule is attracted to its image . When N 
is large enough this attraction becomes more important 
than the loss of entropy due to adsorption. Thus, at a 
given pH>pI the charge of the complex changes its sign 
as a function of growing N. This phenomenon is obvious 
on both phase diagrams of Fig. 1, when at a given pH 
we cross the thick line while increasing N . 

In a more quantitative language, one can say that im- 
age forces lead to the negative correlation chemical po- 
tential of a new PE molecule on the surface of already 
neutralized protein, fJ, s (q, 14). If fj, s (q,lA) < fJ,b{N), 
where Hb(N) is the (negative) chemical potential of a PE 
molecule in the bulk of solution, the negative bare protein 
becomes overcharged by adsorbed positive PE molecules 
so that the whole complex becomes positive. This can 
happen at pH<14 as well. 

Charge inversion by a PE is used in the gene delivery in 
order to invert the negative charge of DNA by a positive 
PE. This facilitates DNA penetration through a nega- 
tive cell membrane || . The change of the sign of DNA 
with increasing N was recorded by the sign change of 
electrophoretic mobility [0. (See also the recent review 
paper on physics of charge inversion j?j and references 
therein). 

Now we can explain the origin of the difference between 
the scenarios A and B. Let us discuss what happens when 
pH increases from pi to 14 for a protein with a given q. 
Consider the surface chemical potential of a PE molecule 
adsorbed to a neutral protein-PE complex, fi s (q, pH), or, 
in other words, free energy of adsorption of an additional 
PE molecule by the neutralized complex. We already em- 
phasized that at pH 14 the chemical potential (i s {q, 14) is 
negative. Actually it is natural to assume that fi s (q, pi) 
is negative, too. Indeed, at pH = pi there are q positive 
and q negative charges more or less randomly distributed 
at the surface of the protein globule. A PE molecule can 
be adsorbed at the surface due to spatial fluctuations of 
this charge. Positive monomers of a PE molecule can ap- 
proach preferably negative charges of the surface avoid- 
ing positive ones. Of course, both the intrinsic rigidity of 
PE and the electrostatic rigidity due to repulsion of PE 
charges do not let the PE molecule to avoid all positive 
surface charges. The PE molecule chooses a spacial scale 
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of fluctuations of the surface charge which it can adjust 
to without too much of loss of the elastic energy. As a 
result of this optimization the PE molecule can bind to 
the neutral protein. 

Thus we see that both jJb s {q, pi) and ^ s (q, 14) are neg- 
ative. Now we have to consider two cases: A) n s (q, pi) 
< n s (q 7 14) and B) fi s (q, pi) > ^ s (q, 14). They generate 
the two above mentioned scenarios A and B. 

In the scenario A, when A and /it, (A) grow, at the 
concentration N = N\, where /Z{,(Ai) = /J, s (q, pi), the 
first PE molecule is adsorbed by the protein. There- 
fore, at A > Ni we deal with the isoelectric point of 
the protein-PE complex, pl(q,N). Positive charge of 
adsorbed PE molecules shifts pl(g, A) to pH>pI in or- 
der to add the compensating negative charge to the pro- 
tein itself. The number of adsorbed PE molecules grows 
with A and pushes pl(q, A) to larger and larger values. 
At N — N 2 , when /zt,(A 2 ) = n s {q,lA) isoelectric point 
pl(<7, A) reaches its upper limit. 

Let us switch to the opposite case fi s (q, pi) > ^ s (q, 14), 
which corresponds to the scenario B. Consider what hap- 
pens in this case with the charge of the protein complex 
when A grows (see Fig. 1). There are again two charac- 
teristic concentrations Ai and A 2 , defined by the same 
equations /it,(Ai) = fJ- s (q, pi) and )J,b(N 2 ) = ^ s {q, 14). In 
the scenario B, Ai > N 2 , and therefore, there are three 
different ranges of A: N < N 2 , N 2 < N < Ai and 
A > Ni. 

1. A < N 2 . If A is so small that (J, b (N) < fi s (q,14), 
then there is no charge inversion of the protein by PE 
and, therefore PE does not affect the focusing. The pro- 
tein is positive at pH < pi and it is negative at pH > 
pi. 

2. N 2 < A < Ai. In this range fi s (q, 14) < fj, b (N) < 
/i s (<7, pi). When pH exceeds pi by some finite value, 
the negative protein starts to adsorb PE molecules but 
the protcin-PE complex still remains negative. At pH 
= pl(q,N), where /i s (<7, pH) = /if, (A) the charge of 
protein-PE complex goes through zero and becomes pos- 
itive (charge inversion). At pH > pl(q,N) the protein 
charge stays positive. 

3. A > Ni. For pH > pi in this case, Hb(N) > fj, s (q, 
pH) so that protein adsorbs PE and has positive charge. 
For pH < pi the protein is always positive, too. There 
are no isoelectric points and no focusing. 

We see that there is only a limited range of N 2 < N < 
Ai, where we get the new (/-dependent isoelectric point 
pl(q,N), so that the suggested separation method can 
work. Note that for each direction of electric field one 
of the two isoelectric points is stable, while another is 
unstable. In the former case on the both sides of the 
point proteins move to the isoelectric point, in the latter 
one they move away from it. The standard isoelectric 
focusing can reveal only one of the two focus points. For 
example, if positive proteins with different q start on the 
left side of the device and drift to the right, i. e. in the 



direction of increasing pH (the electric field is directed to 
the right), all proteins focus in pi. If we reverse the elec- 
tric field so that the positive protein-PE complexes drift 
from the right side of the device to the left, different 
proteins with the same pi focus in the new (/-dependent 
isoelectric points, pl(g, A). Thus we can use this depen- 
dence on q to separate proteins with different q, but with 
the same pi. 

There is however a way to see both isoelectric points 
for a fixed direction of electric field. Let us assume that 
in the absence of electric field proteins are uniformly dis- 
tributed in space with pH gradient. When we apply elec- 
tric field, they concentrate in the stable point and escape 
from the unstable isoelectric point. As a result we should 
observe the maximum of the protein concentration in the 
former point and the minimum in the latter one. 

Until now our theory was strictly phcnomcnological. 
Actual calculation of chemical potentials of fi s (q, pi) and 
Hs{q, 14), which helps to choose between the scenarios A 
or B, requires a detail theory, which takes into account 
the distribution of protein aminoacids, the shape of the 
protein globule, flexibility of PE and the protein, and 
their geometrical dimensions. 

Below we give some microscopic estimates of /i s (q, pi) 
and Hs{q, 14) in a toy model, which is not reliable enough 
to choose between scenarios A and B for real proteins and 
PE (although it seems that the scenario B is more likely). 
These estimates, however, help to understand physics of 
fi s (q, pi) and n 8 (q, 14) and to explain the origin of their 
dependence on q, which leads to possibility of protein 
separation. 

Let us consider a toy model of a protein of approx- 
imately 200 aminoacids as a rigid sphere with radius 
R = 2 nm and q = 20 positive and negative charges 
randomly distributed at the surface of the sphere. We 
assume that PE is strongly charged with the linear den- 
sity of the order of rj = e/ls, where e is the proton charge, 
Ib = e 2 jeksT ~ 0.7 nm is the Bjcrrum length, e ~ 80 
is the dielectric constant of water. Assume that PE has 
Z = 4 charges, so that its length (Z — 1)Ib is comparable 
to R. 

We concentrate on the estimates of zt s (<z, pi) and 
(i s (q, 14) and their comparison. Let us first estimate 
/J, s (q, 14). At pH 14 the charge of the globule — q is large, 
so that protein attracts a number of PE molecules. A 
protein with q = 20 is neutralized by five PE molecules 
with the charge Z = 4. The distance between near- 
est adsorbed PE molecules is of the order of 1.5R. Let 
us assume that the screening radius of monovalent salt, 
r s = 1 nm, i.e. it is three times smaller than the distance 
between PE molecules. Then, in the first approximation, 
one can neglect energy of screened repulsion between ad- 
sorbed PE molecules and calculate zt s (q, 14) as the energy 
of attraction of PE molecule to the surface of the sphere 
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(1) 



For q = 20 and Z = 4, Eq. ([[]) gives {i s (q, 14) ~ -Uk B T. 
This estimate disregards the fact that dielectric constant 
of the globule is actually much smaller than 80. This 
leads to the positive image charge of PE inside the pro- 
tein sphere under each PE. Repulsion from the image lifts 
all adsorbed PE molecules above the protein surface and 
diminishes attraction between PE and protein ||. The 
distance of PE from the protein surface can be calculated 
balancing electric field of the protein qe/(R+ d) 2 and of 
the image rj/4d. This gives d ~ 0.4 nm. In Eq. ([!]) for 
H s (q, 14) one should now replace r s by r s — d and also add 
repulsion energy of PE and its image which can be esti- 
mated as (l/8)r]Zelri(r s /d). Together these two changes 
lead to Hs(q, 14) — — 8k B T. At the same time the ne- 
glected above repulsion between adsorbed PE molecules 
adds another positive term. It can be calculated if we 
sum all exponentially screened repulsion energies of near- 
est neighbor PE molecules and then take derivative with 
respect of number of PE molecules. In this procedure 
only nearest neighbors should be taken into account be- 
cause of the exponential decay of the screened potential. 
This calculation gives correction to n s {q, 14) of the order 
of 2k B T, which results in ^ s (g, 14) ~ —6k B T. 

In Eq. ([!]) we clearly see the origin of the g-dependence 
of fJ- s {q, 14): proteins with larger q stronger attract 
PE. This dependence is of course translated in the q- 
dependence of /i s (q, pH) and N\. Thus, all details of the 
function pl(g, N) are strongly g-dependent what makes 
possible to separate proteins with the same pi and dif- 
ferent hidden charges q. 

Let us now estimate fi s (q, pi). First, we show that 
any strongly charged PE is actually quite rigid due to 
the Coulomb repulsion of its charges. Consider the PE 
electrostatic tension force F, which makes a PE straight. 
It depends on the linear density of charge ry, screening 
radius r s and the PE persistence length a. The energy 
of PE of the length L PE is L PE (r] 2 /e) ln(r s /a). Taking 
derivative with respect of Lp B we get 



F 



(2) 



Now we should recall that the density of protein surface 
charge fluctuates creating a bending force for a adsorbed 
PE. We assume that we are dealing with random charges 
with the density of charge of each sign q/AirR 2 or the 
average distance between them A — R(Ai: / q) 1 ! 2 . 

In order to find the energy of adsorption of a PE 
molecule in the PI point we should consider compe- 
tition between the gain of the Coulomb energy which 
a PE molecule enjoys when it bends to reach nega- 
tive charges and the loss of the elastic energy neces- 
sary to do that. An important role is played here 
by the so called Larkin length 0, L, which tells how 



long is the section of PE, which moves aside by the 
distance r s to use a typical fluctuation minimum of 
the electrostatic potential averaged over the stripe of 
length L and the width r s . This fluctuation of the po- 
tential equals — (e/er s )(ir s /A 2 ) 1 / 2 (r s /i). Then a sec- 
tion of PE with the length L gets the Coulomb energy 
— {eq/e){Lr s /A 2 ) 1 / 2 and the elastic energy Fr 2 /L. Thus 
the change of the total energy of PE equals 



E 



e?y [Lr s fl 2 Fr 2 
"7 A + ~L 



-T ■ (3) 



Optimizing the right hand side of Eq. (JsJ) with respect of 
L and using Eq. (§) we S et L - r s (ATj/e) 2/3 \n 2/3 (r s /a). 
For q = 20A, A ~ 0.8i? and for i] = e/l B we get large 
L ~ R. This results in quite small E ~ k B T. The reason 
that E is small is the relatively small number of charges 
of the typical protein q (and large distance between them 
A). 

We emphasize, however, that E is the attraction en- 
ergy of a PE molecule to a typical spot on the sphere. 
Actually, a PE molecule can find at the protein surface 
a spot, with anomalously large negative surface density 
of charge, where attraction is stronger. Mentally mov- 
ing a PE molecule along a protein globule surface we can 
study distribution function of binding energies, which has 
a lower energy tail and terminates at the lowest energy 
which we want to find. The shape of the tail (gaussian or 
exponential) and the lowest energy are strongly model de- 
pendent (see similar problems in semiconductor physics 
p0[). The absolute minimum of energy can be larger 
than the energy £ in a typical spot by the logarithm of 
the number of different independent states of a PE on the 
surface of the protein or by the square root of this loga- 
rithm. This factor can reach 3 or 4 so that |/x s (g,pl)| can 
be comparable to \fJ. s {q, 14) | , but most likely it is smaller 
and the scenario B is realized. 

For this scenario the concentration N2 can be as small 
as 10 mM and the concentration N± can be 10-50 times 
larger. It is clear from our estimates that both Ni and 
N2 depend on q and this opens multiple possibilities of 
separation of proteins with the same PI and different q. 
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